Synthetic partial waves in ultracold atomic collisions 
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Interactions between particles can be strongly altered by their environment. We demonstrate a 
technique for modifying interactions between ultracold atoms by dressing the bare atomic states 
with light, creating an effective interaction of vastly increased range that scatters states of finite 
relative angular momentum at collision energies where only s-wave scattering would normally be 
expected. We collided two optically dressed neutral atomic Bose-Einstein condensates with equal, 
and opposite, momenta and observed that the usual s- wave distribution of scattered atoms was 
altered by the appearance of d- and g-wave contributions. This technique is expected to enable 
quantum simulation of exotic systems, including those predicted to support Majorana fermions. 



Interactions can be generally understood as arising 
from collisions between particles which exchange energy 
and momentum. For an interacting gas at the lowest tem- 
peratures T, the thermal de Broglie wavelength AdB can 
vastly exceed the characteristic range £ of the interparti- 
cle potential; the interactions can then be cast in terms 
of an effective contact potential that couples states with 
zero relative angular momentum (s-wave interactions). 
Only at higher collision energies, when AdB ^ A do the 
details of the potential become important. 

For alkali atomic gases at degenerate temperatures 
(T < 1 /iK), AdB ^ 100 nm is much larger than the 
characteristic range of the interatomic van der Waals po- 
tential pQ, ^vdw = 1-5 nm to 5 nm, and s-wave inter- 
actions dominate. Collisional Fano-Feshbach resonances 
PQ, which are used to modify the strength of atomic in- 
teractions, do not change the characteristic range of the 
potential and only s-wave interactions are usefully mod- 
ified 0. 

Interactions with higher order partial waves are crit- 
ical for creating systems whose elementary excitations 
are Majorana fermions, enigmatic particles that are their 
own anti-particles. Some many-body systems, includ- 
ing 2D p-wave superfluids [3] and ID Fermi gases with 
spin-orbit coupling [4 , have been theoretically shown to 
support Majorana fermions among their elementary ex- 
citations. Ultracold Fermi gases could realize these sys- 
tems using laser-induced p-wave interactions [5, 6 , the 
fermionic analog to the bosonic d-wave interactions we 
demonstrate here. 

We effectively screen the native atomic interaction with 
light, reminiscent of effects familiar in condensed matter 
systems, such as the screening of the Coulomb interac- 
tion between a pair of electrons by the collective response 
of the remaining electrons [7 . Similarly, the laser dress- 
ing modifies collisions between dressed atoms producing 
effective higher order partial waves, even though atoms 
still microscopically collide in the s-wave channel. The 



interaction potential between a pair of particles located 
at positions r and r' can often be described as a function 
of their separation, V (r — r'). The two-body term 

x V-(k3-ki)^(k 2 )^(ki) 

X 5™ (k 3 + k 4 - ki - k 2 ) , (1) 

in the Hamiltonian describing this interaction is ex- 
pressed in terms of the Fourier-transformed potential 
V(Ak), which describes the transition amplitude for a 
collision that changes the wave vector of a particle by 
Ak = k3 — ki. The field operator </>(k) denotes the anni- 
hilation of a particle with momentum ftk = h(k x , k y ,k z \ 
and the Dirac ^-function ensures conservation of momen- 
tum. The bare interaction between ultracold atoms is 
well described by a contact potential V(r— r') w gS^ (r— 
r'), with Fourier transform V'(Ak) = g = 47rft 2 a s /m, 
where a s is the scattering length and m is the atomic 
mass. As V'(Ak) is independent of Ak, low-energy scat- 
tering is isotropic. The modification of interactions by an 
environment (screening) can be characterized in terms of 
a response function [7] x, producing the screened inter- 
action V' = x(ki, k 2 , k 3 , k 4 )V"(k 3 - k x ). 

To modify the effective atomic interactions we il- 
luminated a Bose-Einstein condensate (BEC) of 87 Rb 
atoms with a two-photon Raman field, producing dressed 
atoms — spin and momentum superpositions. Two A = 
790.1 nm Raman laser beams counterpropagating along 
±e x coupled the three Zeeman states comprising 87 Rb's 
F = 1 ground state manifold with strength HQr, 
(Fig. 1A). The dispersion relations of laser-dressed atoms 
(Fig. IB) are altered from those of free particles, and con- 
sist of three distinct bands labeled by n = { — 1,0, +1} 
with energies i?o,±i(k). These altered single-particle dis- 
persions lead to synthetic magnetic [8] and electric fields 
[9], and spin-orbit coupling [10]. The current work re- 
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FIG. 1. Creation and detection of higher order partial waves. 
(A) Two Raman laser beams couple the three hyperfine Zee- 
man states comprising the F — 1 manifold of the 87 Rb 
5Si/2 electronic ground state. (B) Energy-momentum disper- 
sion relations for Raman-dressed atoms with Raman lasers 
counter-propagating along ±e x . (C) Schematic represen- 
tation of two Raman-dressed BECs with atomic momenta 
hk x = ±2hkL following their collision, indicating the scatter- 
ing halo as imaged after TOF. A magnetic bias field, Boe y , 
generates the linear (lj z ) and quadratic (e) Zeeman shifts. 



veals how Raman-dressing modifies the interaction be- 
tween atoms. For atoms in the n = — 1 (ground) band, 
we demonstrate higher order effective d- and g-wave scat- 
tering between ultracold atoms at collision energies or- 
ders of magnitude below those traditionally required to 
depart from the s- wave regime. Odd partial waves are 
absent when elastically scattered atoms lie in the same 
band and are hence indistinguishable, whereas inelastic 
(band-changing) transitions can have contributions from 
all partial waves. 

We first consider the elastic interactions between 
dressed atoms in the n = — 1 band when the relevant en- 
ergies (kinetic, interaction, and thermal) are less than the 
« HQr energy gap between bands. The dressed states are 
related to the bare atomic spin states by a momentum- 
dependent unitary transformation, U (k) (a 3 x 3 matrix 
which depends on k x , not k y or k z ), that links the spin 
character of an atom to its momentum. To a very good 
approximation, interactions between 87 Rb atoms in the 
F = 1 manifold are spin-independent [11] and cannot 
change the spin, hence the scattering amplitude from one 
pair of momentum states to another depends on the over- 
lap of the spin wavefunctions associated with the initial 



and final momentum states. While microscopically the 
atoms still collide in the s-wave channel, the result of this 
momentum-dependent overlap is that atoms experience 
an interaction altered [12] by a response function 

x = (n=-l\U(k i )W(k 2 )\n=-l) 

x (n = -l|t/(k 3 )t/ t (k 1 )|n=-l), (2) 

leading to scattering described by effective higher order 
partial waves. The modified scattering can be described 
by an effective real space potential (determined by the 
momentum-dependence of \) that is finite-ranged along 
the direction of the Raman beams, as discussed in detail 
in [12]. 

We investigated interactions between Raman-dressed 
atoms by observing the scattering halo formed by the col- 
lision of two BECs in the lowest energy Raman-dressed 
band (Fig. 1C). The experiments started with a nearly 
pure 87 Rb BEC of around 5 x 10 5 atoms in a crossed 
optical-dipole trap with frequencies, {uj x , uo y , uj z } /2ir — 
{13,45,90} Hz. That BEC was then prepared in the 
n = — 1 band, with coupling strength Qr, and split 
into two (initially spatially overlapping) condensates with 
momenta hk x = ±2tt^ per atom [12 . After this 
preparation the dipole trap was immediately turned off 
(t ff < 1 fjs) allowing the dressed condensates to col- 
lide in the absence of any confining potential. The scat- 
tering halos formed in « 2 ms during which time Qr 
was held constant and the colliding BECs separated and 
expanded. After this initial 2 ms stage, we transferred 
the atoms from the ground dressed state into the bare 
\F = l,rap = +1) spin state [9] [12], and imaged the 
atomic distribution after a 36.2 ms time-of- flight (TOF). 

Figure 2 shows representative scattering halos from 
colliding BECs, with (Fig. 2, A-D) and without (Fig. 2E) 
Raman dressing. Both the absorption images of the 
scattering halos (which record the atomic column den- 
sity in the e^-e^ plane) and the radial density distri- 
butions n(k p ,k x ), reconstructed using the inverse Abel 
transform [12] . are depicted. A striking difference be- 
tween the images in Fig. 2, A-D and that of Fig. 2E 
is the shape of the scattering halo. The halo repre- 
sents a constant energy surface on which the scattered 
particles lie as a consequence of conservation of energy 
and momentum in the collision process. For atoms 
with an isotropic dispersion relation, the scattering halo 
is approximately spherical in TOF (Fig. 2E). However, 
the Raman lasers break the spherical symmetry of the 
system by modifying the dispersion relation along e x , 
ELi(k) = S-i(k x ) + h 2 (ky + k 2 z )/2m, resulting in an ap- 
proximately ellipsoidal equipotential surface. 

As VLr increases, £Li(k) tends to the isotropic 
ft 2 k 2 /2m free-particle dispersion, and the scattering halo 
becomes spherical (seen in the progression in Fig. 2, A- 
D). A quantitative analysis of the scattering halos' shape 
is considered in [12]. We restricted our measurements to 
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FIG. 2. Scattering halos produced by colliding Raman-dressed BECs. 2D momentum distributions, projections of the scattering 
halo onto the e x -e y plane, are shown along with the radial momentum distributions, n(k p ,k x ), reconstructed using the inverse 
Abel transform, where k p = (k 2 , + k 2 ) 1 ^ 2 . Darker colors correspond to regions of higher atomic density. (A-D) The halos 
produced by colliding dressed condensates are nearly ellipsoidal and have a non-isotropic angular density distribution. For 
reference, (E) depicts the s-wave scattering halo from a pair of colliding BECs without Raman coupling in the \F = l,mj7 = +1) 
spin state. Each image represents the average of around 50 experimental shots 



HQr > 4 El; for smaller coupling strengths the lowest 
band is significantly anharmonic. The deviations from 
sphericity result from a single-particle effect - the altered 
dispersion relations - and are not a signature of beyond 
s- wave scattering. 

Along with the modification of scattering halo shape, 
the distribution of atoms on the scattering surface is also 
changed. Two main effects determine the distribution of 
atoms: (i) the anisotropic final density of states (DOS) 

(k), and (ii) effective higher-order partial waves in the 
binary collision process. The first effect, resulting from 
the shape of the scattering halo as determined by i£_i(k), 
would occur in an identical scattering experiment with 
spinless bosons subject to the same energy-momentum 
dispersion relation. The second effect, described by the 
response function x °f Eq. |2j is a consequence of the 
momentum-dependent spin structure of Raman-dressed 
atoms leading to effective beyond s- wave interactions, 
and is the central result of this paper. 

These independent effects may be quantified 
with Fermi's golden rule (FGR). For two collid- 
ing atoms the scattering rate from the initial state 
\i) = / i 1 (2/ci / e :E )0 / i 1 (— 2kL^ x ) |vac) into the fi- 
nal state \f) = (^(k/^'i^-k/O | vac) is T(k / ) = 
27r\(f\H[ n t\i}\ 2 pf(kf)/h, dependent on both the final 
DOS pf and the transition matrix element (f\H- int \i) 
describing the screened potential V' . Here |vac) is the 
many-body vacuum, and </4(k) describes the annihi- 
lation of a Raman-dressed atom with wave vector k. 
The dependence of the matrix element on the response 
function x gives rise to effective higher order partial 
wave scattering [12 . The contributions of the DOS 
and the matrix element to the anisotropic nature of the 
scattering rate r(kj) are of comparable magnitude. 



Figure 3A shows the peak signal as a function of polar 
angle extracted from the inverse Abel transform for sev- 
eral different Qr, normalized to the anisotropic DOS [T2] . 
This normalization removes all anisotropics except those 
resulting from beyond s-wave scattering, which we quan- 
tify by fitting the data with the partial-wave expansion 
m, P(f) = A\ E J=0 , 2) 4 (2* + l)(ex P 2% - l)P,(cos0)| 2 / 
/c 2 , where the fit parameter, 7^, is the phase shift of the 
/ th partial wave, P/(cos#) is the Legendre polynomial of 
order /, and A is an overall scaling factor (e.g., depen- 
dent on the initial condensate densities and their inter- 
action time). The peak signal p(0) is proportional to 
T(kj). The fit was restricted to the angular range tt/ 
5 < < 47r/5 which excluded the regions close to the 
unscattered clouds. These regions showed an increase in 
atom number beyond that expected, even for pure s- wave 
scattering between bare spin states, possibly because of 
multiple scattering events or Bose stimulation. 

Figure 3B displays the d- and g-wave amplitudes 
(normalized to the s-wave amplitude) as a function 
of Qr, obtained from fits as in Fig. 3 A, where the 
fractional amplitude of the / th partial wave is (2Z + 
l)sign(7/ ^) |exp2z7# - 1| / |exp2^ - 1|- At Ml R = 
A.IEl, we find a fractional d-wave amplitude of 13% for 
collision temperature T co \ = pv 2 el /2kB ~ 0.7 /iK, where 
p is the reduced mass of the colliding particles and v re \ 
their relative velocity. By comparison, two previous ex- 
periments [HI [15], which studied the scattering of 87 Rb in 
a bare spin state, did not find comparable d-wave mag- 
nitudes until T co i > 100 /iK. The solid lines show the 
theoretical prediction for the fractional partial wave am- 
plitudes, calculated from a partial wave expansion of the 
matrix element (f\Hi nt \i). Figure 3B highlights that the 
engineered effective interaction potentials are not only 
long-ranged but also tunable; by changing we con- 
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FIG. 3. Beyond s-wave scattering. (A) Peak signal of 
the radial momentum distribution (after normalizing to the 
anisotropic density of states) as a function of (circles) and 
fits to a partial wave expansion containing s-, d- and g-wave 
contributions (solid lines). (B) d-wave amplitude (red cir- 
cles) and g-wave amplitude (green circles) as a fraction of 
the s-wave amplitude together with the respective theoret- 
ical values (solid lines). Inset to (B): theoretical fractional 
partial wave amplitudes plotted over a larger range of Qr. 
For HQr/El > (Sk x /kL) 2 (dashed vertical line), the response 
function \ ~ * 1 an d the interactions return to pure s-wave. 



trolled both the magnitude and sign of the d-wave phase 
shift. 

We now consider collisions between atoms in the high- 
est energy n = 1 Raman-dressed band. Although spin- 
changing collisions of bare atoms hardly occur in 87 Rb, 
strong band-changing transitions of dressed atoms are 
in general permitted [I6j [17]. To study this effect we 
investigated the (energetically allowed) collision-induced 
transitions of atoms to lower bands. We considered colli- 
sions in a single BEC (so atoms collided at extremely low 
relative momenta) and observed the spatial distribution 
of the collision products after TOF. 

We prepared a BEC in the highest energy dressed state 
with zero momentum [12 . Figure 4B shows atoms in 
the n = (intermediate energy) band which have de- 
cayed from the excited band by the processes indicated in 
Fig. 4A. Scattered atoms exhibited a node in the k x = 
plane, a dramatic manifestation of beyond s- wave scat- 
tering. The inner halo (I) is due to a process where one 
atom in a colliding pair undergoes an inelastic transition 
to the n = band while the other scattered atom re- 
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FIG. 4. Collision-induced decay from the n — 1 excited 
Raman-dressed state. (A) Decay paths from the n — 1 ex- 
cited state to the n — state (here MIr = 13 El)- (B) 
Momentum distribution of decay products in the n = state. 
After loading into the n — 1 band the atoms were held in 
the dipole trap for 1 ms before release, allowing a significant 
fraction to decay to the lower energy bands. An adiabatic 
deloading technique transferred n — dressed atoms into the 
vriF — spin state [12], and the atoms were subsequently 
imaged after 7 ms TOF. 



mains in the n = 1 band. As the scattering products are 
distinguished by their band index, the spatial wave func- 
tion need not be symmetric with respect to exchange, 
and the scattered wave function can in principle have 
contributions from all partial waves. However, only ef- 
fective p-wave and higher odd partial wave components 
are permitted due to the underlying symmetry of the 
Hamiltonian. 

Our technique for modifying atomic interactions has 
significant ramifications for systems of ultracold atoms, 
bosons and fermions alike. Possible new effects induced 
by synthetic higher-order partial wave interactions in a 
harmonically trapped BEC could include roton features 
in the excitation spectrum [16] , crystalline ordering [18] , 
supersolidity [19 , or even the support of 2D solitons [20] . 
In addition, collisions in this experiment are subject to 
a new kind of resonance phenomena. As the collision 
energy approaches the inter-band spacing, virtual excita- 
tions to the higher band can greatly modify the collisional 
properties beyond the effects described here. 
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METHODS AND MATERIALS 

Dressed states 

Consider the operator a (q) which destroys a par- 
ticle with momentum ftq and a = 0,=bl, where a 
denotes mp = 0, ±1. To describe the states cou- 
pled by the Raman dressing it is convenient to intro- 
duce the operators {0_i(k), (/>o(k), 0+i(k)} = {<£>_i(k + 
2kpe x ) 1 ^o(k), 0+i (k — where the quantum 

number hk is sometimes described as the "quasimomen- 
tum," in analogy to the appropriate quantum number for 
dressed (Bloch) states in a lattice potential. From here 
on and throughout the paper we drop the tilde notation 
and refer to hk as the "momentum" of a dressed parti- 
cle. Our experimental system (Fig. 1) is described by the 
Hamiltonian 

/d 3 k - 
—-3 4(k)i^'(k)<Mk) + JTint, (SI) 

a,a f 

where H a:Cr ' is an element of the 3x3 matrix H = 
Hi(k x ) + [H 2 (ky + k 2 z )/2m] ® i and H int describes bi- 
nary interactions between the particles (discussed be- 
low). The free motion along e x , together with Zeeman 
shifts and atom-light coupling, is described by Hi(k x ), 
1 is the 3x3 identity matrix. (We do not consider 
an external potential as the condensates collided after 
the optical dipole trap had been turned off and the 
scalar light shift induced by the Raman beams was 
made negligible by operating at the magic wavelength 
A = 790.1 nm). A static bias magnetic field Bq e y with 
Bo ~ 406 fiT generated a uz/2tt = 3.25 MHz Zeeman 
shift between the bare spin states mp = 0, ±1, with a 
small he = 0.42 El offset of the mp = state due to the 
quadratic part of the Zeeman shift {E^jh = 3.6 kHz). 
The counterpropagating Raman beams had orthogonal 
linear polarizations and a frequency difference Acjl, (giv- 
ing a detuning S = Aujl — oj z from resonance) cou- 
pling the bare states \rriF = — 1, k x + 2k l), \ttif = 0, k x ), 
\rriF — +1, k x — 2ki}j with strength Qr. In this bare ba- 
sis and in the rotating frame Hi(k x )/h has the form [21] 

^(k x + 2k L ) 2 -S Sl R /2 \ 

n R /2 Ji-(k x -2k L ) 2 +5j 

For the experiments described in this paper (5 = ex- 
cept for a final deloading step (discussed below) after 
the condensates had collided (the deloading procedure 
mapped dressed eigenstates, which are spin-momentum 
superpositions, back to a single bare spin-momentum 
state). For a given value of k Xl Hi can be diagonal- 
ized by a 3 x 3 unitary matrix, ?7(k), and the Raman- 
dressed eigenstates are related to the bare states by the 
momentum-dependent unitary transformation: </4(k) = 
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^Zcr ^n,cr(k)(/>cr(k), where n is the band-index of the 
dressed states. 



System preparation 

Our experiments began with nearly pure 87 Rb BECs 
of around 5 x 10 5 atoms in the \F = l,rap = —1) state 
confined in a crossed optical dipole trap. The trap con- 
sisted of a pair of 1064 nm laser beams propagating along 
e x (1/e 2 radii of w y ~ 120 //m and w z ~ 50 jam) and 
e y (1/e 2 radii of w x « w z « 160 /am). The trap fre- 
quencies at the end of an evaporation sequence were 
{<jj x ,<jjy,u} z }/2ir = {13,45,90} Hz. 

To study collisions beyond s-wave we prepared col- 
liding Bose-Einstein condensates in the lowest energy 
Raman-dressed band using the following procedure: (i) 
starting with an \F = 1, = —1) BEC at rest, a radio- 
frequency magnetic field, B T f(i)e x , was ramped on in 
100 ms with an initial detuning h5 r f = —26 El (where 
S r f = uj r f — wz, and the detuning 5 r f was controlled by 
changing Bq) and coupling strength Ml T f = 3 El; (ii) 
the detuning was adiabatically ramped to 5 r f ~ in 
25 ms, and then the radio-frequency coupling strength 
fi r f was decreased to zero in 12 ms. In the last millisec- 
ond of the ramp of f} r f the dressed state to which we 
connect adiabatically is a virtually equal superposition 
of vfip = d=l; (hi) During this final millisecond Qr in- 
creased from zero to its final value, creating two spatially 
overlapping condensates in the ground dressed band with 
momenta ±2^/c^. The turn-on of the Raman beams was 
adiabatic with respect to inter-band excitation but fast 
with respect to the trap's ~ 20 ms quarter-period along 
e x . This ensured the atoms could not change their mo- 
mentum as was increased [2T] . 

To observe s-wave scattering we prepared the system 
in the absence of any Raman beams. The BEC, initially 
formed in the \F = l^mp = —1) state, was transferred 
to the \F = l,ra.p = +1) state by radio- frequency adia- 
batic rapid passage. Thus the atoms were in the same 
spin state as for the time-of-flight (TOF) stage for the 
higher partial wave data (see Deloading below), ensuring 
that any small distortions of the scattering halo in TOF 
expansion due to stray magnetic field gradients were the 
same for the s-wave data and the higher order partial 
wave data. Two temporally square pulses of a ID optical 
lattice along e x (formed with A = 792.1 nm light) were 
used to split the BEC, initially at rest, into ±2ft/c£ orders 
[22] . The dipole trap was switched off (t ff < 1 /^s) at the 
end of the second optical lattice pulse, and the result- 
ing atomic distribution absorption-imaged after 36.2 ms 
TOF. 

To prepare a BEC in the highest energy dressed state 
we first loaded the BEC into the lowest energy Raman 
dressed state with zero momentum [9j [21] . The BEC was 
then excited to the highest energy dressed state with zero 



momentum by amplitude modulating the Raman beams. 
For a coupling strength Qr = 13 El, the amplitude of 
the Raman beams was modulated by 17% (half peak- 
to-peak amplitude) for 1 ms. During this 1 ms the fre- 
quency of the modulation was swept across the resonant 
energy difference AE between the ground and excited 
dressed states at zero momentum, from AE/h — 3 kHz to 
AE/h + 3 kHz, allowing near complete population trans- 
fer to the highest energy dressed state by adiabatic rapid 
passage (AE = 18.9 E L = h x 69.4 kHz for Q R = 13 El). 
The 1 ms sweep time represented a compromise between 
being slow enough for adiabatic transfer while being fast 
enough to avoid significant decay from the excited state 
during the transfer. 



Inverse Abel transform 

We used the inverse Abel transform [23] to reconstruct 
the radial density distribution, n(p, z), from the column 
densities n(x,y) recorded in our absorption images. The 
inverse Abel transform is appropriate for axially sym- 
metric distributions, which is the case for the scattering 
products after TOF. The small central cloud of atoms 
in the images (Fig. 2) is due to a residual population of 
k = atoms remaining from the loading process into the 
ground Raman-dressed band. The density of the k = 
clouds is small enough that their effect on the scattering 
halos is negligible. 



Deloading 

Raman-dressed atoms are in a superposition of spin 
and momentum states. When the Raman beams 
were turned off diabatically, the Raman-dressed atoms 
were projected onto a superposition of bare spin states 
\m F = -l,k x +2k L ), \m F = 0,k x ), \m F = +l,k x - 2k L ), 
somewhat complicating the atomic distribution in TOF. 
In order to image a single, well-formed scattering halo, 
we adiabatically mapped atoms in the n G { — 1,0, +1} 
band to the \ttif = —n, k x + 2/cz,n) bare state. This was 
achieved by ramping the detuning of the Raman beams 
from HS = to HS = 50 El (for which the Raman-dressed 
states coincide nearly perfectly with pure spin states). 
The 2 ms long ramp was started after 2 ms of TOF, by 
which time the colliding condensates had spatially sep- 
arated. After the ramp Qr was turned off in 100 /is, 
leaving the atoms to continue their time-of-flight in the 
rriF = +1 state. In analyzing the data we neglect the 
change in group velocity during the few milliseconds TOF 
during which the atoms are dressed. This has a negligible 
effect. 
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Effective higher order partial waves 

We now consider the term describing binary interac- 
tions between particles, #mt- The four-field interaction 
term for spin-independent (a very good approximation 
for 87 Rb atoms [IT] but not necessarily for other atomic 
species) contact interactions of the form V(r — r') = 
gS(r — r') is 

Hint =f Jd 3 v $U ( r )^ (r)^ a (r) (S2) 



fiQ R =4.2 E L 



/ 



d 3 ki d 3 k 2 d 3 k 3 ri 3 k 4 - . . 

(27T)3 (27T)3 (27T)3 (27r) 3 V/ ( k l' k 2' k 3,k 4 j, 



with 



V-(k 1 ,k 2) k 3) k4)=| E ^ 1 (k 4 )^ 2 (k 3 )4 1 (k 2 )4 2 (ki) 

x ^ 3 )(k 3 + k 4 -ki-k 2 ). (S3) 

Eq. S4 shows that the matrix element for a collision be- 
tween the bare states is independent of the momenta of 
the colliding particles. 

Transforming to the dressed state basis the interaction 
Hamiltonian becomes 



V 



9 



E 



l(k 4 )^L(k 3 )^ 1 (k 2 )^(k 1 ) 



n\ ,n 2 ,n i ,7i2 



E c/ ™' 1 ^ 1 ( k 4)^ 1) „ 1 (k2)f/„ 2 , CT2 (k 3 )^ 2) „ 2 (k 1 ) 



5 (3) (ka 



k 2 ). 



(S4) 



The key point is that the transition amplitude for a pro- 
cess in which dressed particles with initial wave vectors 
ki and k 2 are scattered to final wave vectors k 3 and 
k 4 is now dependent on k x ,k 2 ,k 3 ,k 4 through the ele- 
ments of the unitary matrices. For example if the initial 
and final scattering states are constrained to lie within 
the lowest band, we have \i) = /t _ 1 (k 1 )0 /t _ 1 (k 2 )|vac), 
\f) = ^ / Li(k3)0 / Li(k4) | vac), and the matrix element be- 
tween these states is 

(f\H in t\i) = |[x(ki,k 2 ,k 3 ,k 4 ) + x(k 2 , ki, k 3 , k 4 ) 
+ x(ki, k 2 , k 4 , k 3 ) + x(k 2 , ki, k 4 , k 3 )] 
x^ 3 )(k 3 + k 4 -k!-k 2 ), (S5) 

where 

X(k 1; k 2 ,k3,k4)= E ^-i,a 1 (k 4 )C/] 1 ,_ 1 (k 2 ) 

xf/. lil! (k 3 )(/t rl (k 1 ). (S6) 

Moreover, it is apparent from Eq. [S4] that the matrix 
elements for band-changing collision events are in gen- 
eral non-zero, and also dependent on the initial and final 




jt/4 jt/2 3jc/4 n 
Polar Angle, 9 

hQ R = 1 1 .8 E L 




it/4 Jt/2 3jt/4 Jt 
Polar Angle, 9 

FIG. SI. Contributions of density of states (solid green line) 
and higher order partial wave scattering (black circles) to the 
dependence of the measured scattered atomic distribution 
(red triangles). The higher order partial wave scattering sig- 
nal is given by the ratio of the red triangles to the solid green 
line (A) hQ R = 4.2 E L . (B) Htt R = 11.8 E L . 



momenta of the colliding particles. In the lowest energy 
Raman-dressed band such band-changing collisions are 
energetically forbidden provided the collision energy is 
small compared to the gap to the excited bands (as is 
the case in this paper) . For particles in the excited bands 
collision-induced decay to lower energy bands is possible 
as discussed in the main body of the paper and shown in 
Fig. 4. 

The total elastic scattering cross section is given by 
the sum of cross sections of each partial wave contribu- 
tion, however the total collision cross section does not in- 
crease when the atomic interactions are modified by the 
laser dressing. This is because rather than adding d-wave 
and g-wave contributions on top of the s- wave scattering 
component, the interactions are modified such that part 
of the s-wave contribution is renormalized to d- and g- 
wave contributions. The total collision cross section is in 
fact generally smaller than the bare state collision cross 
section, due to the partial (spin-dependent) overlap be- 
tween the initial and final scattering states. The modified 
density of states in the presence of laser dressing also re- 
duces the number of final states for a particle to scatter 
into. 
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Density of states 

The density of states on the scattering halo is purely 
geometric in origin and depends on the coupling strength 
Qr, which determines the shape of the scattering halo 
through the modified dispersion relation E(k). Since 
Qr is measured (see below) for each scattering halo, one 
can calculate the form of the density of states and re- 
move its effect from the measured atomic signal, revealing 
the underlying anisotropic atomic distribution due to the 
modification of the interaction potential by the response 
function x(ki, k2, k3, k^). The uncertainty introduced 
by dividing out the density of states is set by the uncer- 
tainty in Qr, from which the bandstructure and hence 
the density of states are calculated. The Raman coupling 
strength was measured using a Raman pulsing technique 
[21]. The uncertainty (standard deviation) in Qr from 
this measurement was ±5%, translating to a compara- 
ble or smaller uncertainty in the extracted partial wave 
amplitude (determined by carrying out the partial wave 
fitting procedure on data which had been normalized to 
the density of states calculated for Qr ± 5 %). The un- 
certainties in the extracted partial-wave amplitudes of 
Fig. 3B represent the combined uncertainties from the 
partial wave fits to the data (which dominate) and the 
uncertainty in the density of states. 

Fermi's golden rule (FGR) states that the transition 
rate from an initial state |i) to a final state |/) due to 
a coupling H int is T if = 27r\(f\H int \i}\ 2 5 (E t - E f ) /h. 
If one presumes there is some parameter, i.e. the wave 
vector k, over which the final states are uniformly dis- 
tributed then the total transition rate can be written 



Peak Locations 




7i/4 7i/2 371/4 
Polar Angle, 



FIG. S2. Experimentally determined location (circles) of the 
scattering peaks, k — (k 2 p + k^) 1 ^ 2 in Fig. 2, A-D (see main 
text) as a function of polar angle together with theory 
(solid lines). The theory curves have no adjustable parame- 
ters and neglect mean- field-driven expansion. Black - 4.2 El; 
magenta - 5.3 El; green - 9.9 El; red - 12.8 El; blue - s-wave. 
The s-wave halos displayed small deviations from a perfectly 
spherical form, as expected for an initially anisotropic con- 
densate for which mean-field interactions lead to a detectable 
distortion of the halo at the few percent level [26] . 



Shape of the scattering halos 

As Qr increases, E-i(k) tends to the isotropic 
ft 2 k 2 /2ra free-particle dispersion, and the scattering halo 
becomes spherical, as can be seen in the progression in 



2tt 



d 3 k 

(2^) 3 



|(/|^ int |z)| 2 J(^-^(k)). (S7) 



Fig. 2, A-D (main text). Figure S2 quantifies the shape of 
the scattering halos, by plotting the positions of the peak 
in n(k pi k x ) as a function of the polar angle 0. The peak 
locations are in agreement with those expected from nu- 
merical calculations of £"_i(k) using the measured values 
of Q R . 



Considering a function K(#) that gives the momentum 
as a function of on the energy surface -E'(K), the total 
scattering rate can be written 



me)] 




[EM 

mow 



sin<9d<9 



(^int|K)| 2 

V,£(K) 



where K' is the derivative of K with respect to 6. The 
measured peak atomic signal at a given angle, p(9), is 
proportional to T, and hence proportional to the quantity 
in square brackets. Figure [Si] displays the peak atomic 
signal as a function of polar angle (red triangles), as di- 
rectly measured from the radial momentum distribution, 
as well as the calculated contribution from the density of 
states (solid green line). The ratio of these is the con- 
tribution attributed to beyond s-wave scattering (black 
circles). It is this latter signal which is shown in Fig. 3 
A. 



Lifetime in the n = +1 excited Raman band 

We observed a r = 12(1) ms 1/e lifetime for decay 
from the excited state, in reasonable agreement with an 
approximate calculation of the initial decay rate due to 
binary collisions in our BEC based on FGR [24], pre- 
dicting a lifetime r « 9 ms. These excited bands are 
of intrinsic interest, for example in schemes considering 
oppositely-charged particles in synthetic gauge fields or 
proposals for generating non-abelian gauge fields for ul- 
tracold atoms, which rely on two excited dressed states 
(see [25] and references therein). 

In order to measure the lifetime in the n = +1 band 
we excited a BEC in the \n = — l,k = 0) state to the 
\n = +l,k = 0) state as described above. The strength 
of the Raman coupling for this lifetime experiment was 
Qr = 7 El. Figure [S3] shows the number of atoms in 
the excited state as a function of hold time in the opti- 
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FIG. S3. Number of atoms in the excited Raman- dressed 
state as a function of hold time in the optical dipole trap 
(circles). The solid red line is a fit to a double exponential 
function. 



cal dipole trap. For simplicity a double exponential fit of 
the form N n= +i(t) = A\ exp(— t/r\) + A 2 exp(— tjr^) was 
used to the extract the initial decay time of t\ = 12 ms. 
We attribute the slower decay, r 2 = 135 ms, to the de- 
pletion of excited state BEC and subsequent decay of 
the lower density thermal atoms. Note that in the limit 
of very large Raman coupling, HQr ^> El, (which was 
outside our range of available laser powers) the matrix 
element for these band-changing collisions tends to zero, 
and long lifetimes in the excited state should be possi- 
ble. It would also be possible to suppress collisional losses 
from the excited state at lower Qr by modifying the final 
density of states, for example, by the use of an optical 
lattice to open band gaps at the relevant energies. 



Real space representation of the effective interaction 
potential 

That the response function x depends on the initial 
and final momenta of the colliding particles implies that 
an effective real space potential describing the screened 
collisions can no longer be represented by a Dirac 5- 
function but has some finite range. Furthermore, the 
effective interaction potential is long-ranged only along 
the direction specified by the Raman beams (which can 
in general be different from the collision axis). 

In order to qualitatively understand the size of this 
range it is useful to consider two limiting situations. 
Firstly, atoms scattering with a small momentum dif- 
ference, h(ks^ x — k^ x ) <C \/2mMTR, have nearly the 
same bare spin composition. Secondly, particles with 
a large momentum, hk x ^> \/2mhQ J R 1 are unaltered 
by the Raman coupling. In either case, % « 1 and 
the atoms collide as they did absent Raman coupling 
(17). For intermediate momenta Hk x « ^2mhVtR there 
is significant dressing-altered scattering, giving a range 
£r ~ ft/ \/2mhQ,R = ^/EL/M},Rk L , to the effective inter- 
action potential, with £r & 50 nm ^> £ v dw for our exper- 



iment. This expression is valid in the limit HQr > 4i^, 
required for restricting atoms to the lowest band for our 
collision energies. Whereas slowly moving atoms experi- 
ence s- wave, scattering higher partial waves are possible 
at slightly larger velocities. This is because the range of 
the effective potential is much larger than that of the real 
potential, and hence the velocity for which higher partial 
waves become important is proportionally smaller. 

To gain further insight into the nature of the real 
space representation of the effective interaction poten- 
tial a quantitative analysis is necessary. An analyti- 
cal approximation to the real space screened potential 
can be found for he = — 4£x. While this is not the 
case for the experiments described in this paper (where 
he = 0.42£x), this limiting case is useful for pedagogical 
reasons. The form of the real space potential for general 
e can be solved numerically and has a similar form to the 
analytical solution we derive below. 

In the following discussion we consider the component 
of an atom's momentum along the direction of the Raman 
beams, dropping the subscript x, that is k x = k. For 
he = — 4El, Hamiltonian |S2| is 



H 



k 2 l 



(4k-5)F z + QF X , 



(S8) 



where F = (F Xl F yi F z ) is the total angular momentum 
vector, Q = Q/y/2, and we have used dimensionless vari- 
ables Q —> Q/El, k — >• k/kL- Owing to the simple form, 
the problem can be solved using angular momentum al- 
gebra and the screened interaction in the lowest dressed 
state is exactly 



1 + 16^1^3/n 2 




16/cf/O 2 



1/2 



16k 2 /& 



1/2 



1 + Wk 2 k 4 /n 2 



1/2 



l + 16fe|/fi 2 l + 16kj/n 2 



1/2 



As is often the case with screening problems, this effective 
interaction potential is not simply related to a coordinate 
space potential Vflri — r 2 \)- In the following sections, 
we will see when such a potential exists and what it is. 

Perhaps the most unusual aspect of this potential is 
that while momentum is still conserved, Galilean invari- 
ance is not. This breakdown of Galilean invariance is 
common in screening problems where there is a preferred 
reference frame. In this case, rapidly moving atoms are 
Doppler shifted from Raman resonance and the effects of 
the laser coupling vanish. 



Galilean invariant form 

In our experiment, the dressed atoms are slowly mov- 
ing, so we expect a Galilean invariant system. By ex- 
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panding % to second order in the center of mass momen- turn fc C om> we find 
I 



X(ki,k 2 ,k 3 ,k 4 ) 



l + 4(k 1 -k 2 )(k 3 -k 4 )/n 2 



i + 4(fc! - fc 2 ) 2 /n 2 



1/2 



i + 4(k 3 - k A y/& 



1/2 



Q 2 



(S9) 



Because this is only a function of the differences k\ — k<z 
and ks — &4, it has recovered Galilean invariance, but 
only when (k com EL/kL^) 2 (expressed in terms of exper- 
imental variables) is small. In our collision experiment, 

^com — 0- 

Still, this dielectric function is not of the form x(ki — 
/C3), so the resulting potential is not just a function of a 
single distance. This potential has a fascinating non-local 
structure in r-space which goes beyond the current work. 
If we introduce new difference variables, Skis = ki — k% 
and if 13 = ki + k$ (and likewise for k<i and k±) and 
assume that the sum of the initial and final momenta 
if 13 is small, we arrive at the expansion 



1 + 4(5*13/0)' 

K 13 if 24 



O 



n 2 



o 



o 



K 2 

Q 2 



K 2 

11 13 

Q 2 



(S10) 



expressing the interaction in /c-space. This gives rise to 
the screened potential, 



V(Sk) « g 



1 



_l + 8((5fc/0) 2 



(Sll) 



The effective real space potential can be found from the 
Fourier transform of Eq. |S11[ giving 



V(r) 



32 



(S12) 



Recall we made ft — > Q/El, k k/kL, and x — > xkj,. 
The potential is normalized to g = Aith 2 a s /m so it cor- 
rectly represents the gS(x) contact interaction as £1 — >• 00. 
The range of the potential is ~ \/[k(Q/ El)], which 
agrees with the qualitative estimate of « 50 nm for typ- 
ical experimental parameters. 



Effective pseudopotentials 

In the above, we arrived at an effective potential be- 
tween atoms. If instead, we assume all of the momenta 
are small, but keep the second order terms, we can derive 
an effective pseudopotential 



V(x)xsg 



5{x) 



k 2 L n 2 



(S13) 



in which the second-derivative of the delta function leads 
to (i-wave scattering. 



